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1. INTRODUCTION 

Fault diagnosis has become an important topic in industrial applications, which further leads to more 
attention in research community. Model based fault detection, isolation and identification has received much 
attention in the literature [1],[2],[3],[4],[5],[6]. 

Generally, obtaining a proper model of the system is crucial for fault diagnosis design. Presenting a 
proper model by considering multi physics nature and the interaction between continuous and discrete dynamics 
of the system is a challenging task in the model based FD system design for hybrid systems. Bond Graph (BG) 
method may be considered as a convenient tool for modeling and fault diagnosis of dynamic systems. Different 
applications of the BG method can be found in the literature for modeling [7],[8],[9],[10],[11],[12] and fault diag- 
nosis [12],[13],[14],[15],[16] applications. Because of its efficiency in modeling and fault diagnosis of continuous 
dynamic systems, the BG method is developed to represent switching phenomenon in hybrid dynamic systems, 
which leads to numerous applications of the method for modeling and fault diagnosis of hybrid systems in the 
literature. 

Due to the advancement of computer technology in control systems, switching systems have drawn much 
attention in recent years. Switching systems may be assumed as an important class of hybrid systems wherein a 
finite number of dynamical subsystems and a logic rule define the overall dynamic of the system. The dynamic of 
the subsystems and thus, the overall system may be continuous or discrete, linear or nonlinear and so on. In this 
current paper, the linear continuous case of the switching system with average dwell time is considered. 

Fault detection of the switching systems draws much attention in recent years. In [17], fault detection 
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problem for a class of linear switched time varying delay system is formulated as H filtering problem with 
ADT in LMI form. A filter based fault detection approach for continuous time switching systems with ADT is 
considered in [18]. Simultaneous Fault Detection and Control (SFDC) for continuous and discrete time linear 
switched systems with ADT based on dynamic filter is tackled in [19] wherein a fault sensitivity is achieved for a 
given disturbance attenuation level based on H_ /H,. criteria. 

The robustness of an FD system, which is defined as the considering the effects of uncertainties, distur- 
bances and noises in the system, may be regarded as the most important criteria in an FD system. However, it 
is impossible to completely removal the effects of uncertainties, disturbances and noises in an FD system. Thus, 
many studies are devoted to attenuate the effects of disturbances, uncertainties and noises in an FD system. In other 
hand, the attenuation of disturbance effects may lead to reduce the fault sensitivity in an FD system. Therefore, the 
fault sensitivity must also be considered in an FD system design. 

The effects of uncertainties in the model may lead to missed or false alarms in the FD system. Toward this 
end, two approaches are proposed in the literature including active and passive approaches [20]. In active approach, 
the robustness of the FD system is defined in the residual generation stage, while the robustness is considered in 
the residual evaluation stage in the passive approach [20]. 

In [21], the robust fault detection filter design of continuous time switched delay systems is considered. 
Design of observer-based robust power system stabilizers by considering parametric uncertainties is noticed in 
[22]. The parametric uncertainty representation in the BG method in Linear Fractional Transformation (LFT) form 
is presented in [23] and noticed in thereafter studies for robust fault diagnosis [24],[25]. Robust fault diagnosis 
of energetic system with parameter uncertainties is tackled in [24]. The nonlinear modeling, structural analysis, 
residual generation with adaptive thresholds and sensitivity analysis are done using the BG method. The proposed 
method is then used for a boiler system. The modeling and fault diagnosis of a DC motor relying on the BG method 
is given in [16]. The parameters of the BG method are obtained using the real data of the system. The average 
values of the parameters and their standard deviations are assumed as the nominal values and uncertain parts of the 
parameters, respectively. In [25], the residuals and thresholds generation in presence of parameter uncertainties 
in LFT form is considered and is applied for mechatronic systems. Robust fault diagnosis and prognostics of a 
hoisting mechanism based on the BG method considering the LFT form of the parametric uncertainties is addressed 
in [26]. 

Several studies are concentrated on the fault diagnosis of switched system and robust fault diagnosis of 
continuous time systems as well, while there are few works on the robust fault diagnosis of switched systems. 
Incremental BG approach is given for hybrid systems in [27], which is an extension of the proposed method 
in [28] for continuous time dynamic systems. The adaptive thresholds are obtained considering the parametric 
uncertainties in the system, which are dependent on the modes of the hybrid system as well. In [29], a robust 
fault detection and isolation on the basis of pseudo BG model for hybrid systems is presented. The parametric 
uncertainties are given in the LFT form and are assumed 2% of their nominal values at maximum. 

Robust FD systems based on the BG method may be considered as the passive approach, since some 
thresholds are assumed in the FD system to detect or isolate the occurred faults. In this paper, based on the 
BG method, a new active robust FD system is proposed in which the disturbance is attenuated, the effects of 
uncertainties are considered and the fault sensitivity is enhanced as well. The proposed method has the benefits of 
both the BG and the observer method. 

In summary, the main contributions of the paper may be stated as follows: 


e A new robust active fault detection based on the combination of the observer and BG method for linear 
switched system is presented, which simultaneously attenuates disturbance level and enhances fault sensi- 
tivity. 


e Error form of the GARRs is used as the residuals which is based on the output estimation error of the 
observer. 


e The effects of disturbance and parametric uncertainties in the thresholds are considered as a fixed threshold 
and an adaptive threshold, respectively. 


The remainder of the paper is organized as follows. The problem formulation and some lemmas and re- 
marks based on these lemmas are given in section 2. The Research Method of the paper is given in five subsections 
of section 3 as disturbance attention level, enhanced fault sensitivity, simultaneous optimal fault sensitivity and 
disturbance attenuation level, the GARRs in error form based on the BG method and robustness of the FD system 
against parametric uncertainties. The simulation results and discussion for a two-tank system are given in section 
4, follows by conclusion in section 5. 
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2. PROBLEM FORMULATION AND PRELIMINARIES 
Consider the continuous time linear switching system state space representation in the form of 


&(t) = Aux (£) + Buu(t) + Bad(t) + By f(t) (1) 
y(t) = Cizx(t) + Daid(t) + Dp f(t) (2) 

which can be given as Eq.(3) - Eq.(4) based on the BG model of the switching system. 
&(t) = A(a,0)x(t) + B(a,0)u(t) + Ba(a)d(t) + By (a) f(t) (3) 
y(t) = C(a, @)a(t) + Da(a)d(t) + Dz (a) f(t) (4) 


where x € R” is the states, y € RP? is the outputs, u € R is the inputs, d € R? is disturbances and f € R” is faults 
in the switched system. A ,B and C are respectively system, input and output matrices with appropriate dimensions 
denote the system, input and output matrices, respectively. . Ba, By, Da and Dy are the matrices of disturbance 
and fault distribution on the states and on the outputs of the system, respectively. The 7,2 € l = 1,2,..., N index 
is considered to define modes of the switched system. The a = [a, a5 . . . az] is the state of the controlled junctions 
in the BG model, which will define the active mode of the switching system. The number of controlled junctions 
in the BG model is assumed as z. The system matrices are dependent on the parameters of the BG model (R, C, 
L, ...) which represented as 0 in the model. 
The switching Luenberger observer is assumed in the form of Eq.(5). 


a(t) = Ai&(t) + Liy(t) — $(0) (5) 
g(t) = Cist) 
where L; is the observer gain and is dependent on the active mode of the BG model. The observer gain is 


designed for each mode in such a way that some criteria are satisfied. 
The output estimation error, which is defined as Eq.(6), is used in the new form of the GARRs. 


ey(t) = y(t) — g(t) (6) 
The asymptotical stability of the observer is investigated on the error dynamic of the observer as Eq.(7). 
e(t) — (A; x LiCi)e(t) + (Bri = L;D yi) f(t) + (Bai + LiDqai)d(t) = (7) 

Aaie(t) + Bagif(t) + Baaid(t) 


The output estimation error is achieved as Eq.(8). 
€y (t) = C;e(t) + Dyif (t) + Daud(t) (8) 


In this paper, two criteria are considered for the output estimation error of the observer for disturbance 
attenuation level and fault sensitivity, which is presented in Eq.(9) and Eq.(10), respectively. 
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The following definitions and lemmas are used in the current paper. 

Definition. Let N5(t1,¢2) stand for the switching number of ó(t) on the interval [t1, t2) for a given 
switching signal ó(t) and any t2 > tı > to. If Eq.(11) is satisfied for No > 0 and Ta > 0, then 7; is called the 
ADT and Nj is chatter bound [30]. 


Na(t1,t2) € No + (te — t1)/Ta (11) 


Lemma 1: For the continuous time switching system Eq.(1) and Eq.(2), let a > 0, u > 1 and y; > 
0, Vi € I be constant scalars. Suppose there exists positive definite C! function V; : R” — R,ô € l with 
Vs(to) (x (to)) = 0 such that [31]: 


Vi(z;) € —oVi(zi) — yT + yulu Ni €1 (12) 
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Vi (tk) < UV; (zu), V(t, j) Elx li Fi 


then the switching system is Globally Uniformly Asymptotically Stable (GUAS) and will satisfied Eq.(9), which 
is called Ha performance, with index y no greater than max(»;) for any switching signal with ADT as Eq.(13). 


Ta > T} = In(y)/o (13) 


V is called the Lyapunov function. 
Remark 1: In order to use the Lemma 1 in this paper, the u and y are substituted by d(t) and e, (t), 
respectively, which indicate the disturbance attenuation level on the output estimation error of the observer. 
Lemma2: For the continuous time switching system Eq.(1) and Eq.(2), let a > 0, u > 1 and B; > 
0, Vi € lbe constant scalars. Suppose there exists positive definite C! function Vsa : R” > R,ô € l with 
V5 (to) (@(to)) = 0, which is called Lyapunov function, such that [32]: 


Vi (22) € —aVi(zi) + ye yt — Butu Vi El (14) 


Vi (tk) < UV; (zu), V(i, j) Elx li FJ 


then the switching system is GUAS and will satisfied H_ performance as Eq.(10) with index 5 no smaller than 
min(;) for any switching signal with ADT as Eq.(13). 

Remark 2: In order to use the Lemma 2 in this paper, the u and y are substituted by f(t) and e, (t), 
respectively, which indicate the fault sensitivity on the output estimation error of the observer. 

Lemma 3 [33]: For a given m x m symmetric matrix Z € Sm and two matrices U and V of column 
dimension m, there exists matrix X that is unstructured and will satisfy Eq.(15). 


UTXV -V"XTU -Z «0 (15) 
if and only if the following inequalities as Eq.(16) and Eq.(17) with respect to X are satisfied. 


NUZNyu <0 (16) 


NUZNy «0 (17) 
Ny and Ny are arbitrary matrices that their columns are a basis for null spaces of U and V, respectively. 
3. RESEARCH METHOD 


In this section, the LMI formulation for fault detection problem of linear continuous time switched system 
Eq.(1)-Eq.(2) is given. The overall results are given in five subsections. 


3.1. Disturbance attenuation performance in the output estimation error of the switched observer 


The switching Lyapunov function for disturbance attenuation criterion of the observer is considered as 
Eq.(18), which its derivative must be negative definite as Eq.(19). 


V(t) = e(t)? P;e(t) > 0 (18) 


V (t) = e(t)? Pye(t) + e(t)" Pje(t) < 0 (19) 
The switching Lyapunov function derivative may be obtained as Eq.(20) in fault free case. 
V(t) = (e(t) AS; + A(t)” Boras) Pie(t) + e(t) Pi(Aaie(t) + Baaid(t)) (20) 
which can be written as Eq.(21) for simplicity. 
V(t) = e(t) (AT; P; + PiAcu)e(t) + 2e(t)" P; Baaid(t) QD 


According to Remark 1, we need e, (1)? e, (t) for using lemma 1, which can be given as: 


ey (t) e(t) = e(t)! CT C;e(t) + 2e(t)* C7 Da;d(t) + d(t)? Dj, Daid(t) (22) 


IJECE Vol. 8, No. 4, August 2018: 2157 - 2171 


IJECE ISSN: 2088-8708 


Thus, Eq.(12) may be obtained as follows using Lemma 1. 


e(t)" (AT; P; + P;Aci)e(t) + 2e(t)* P; Baa;d(t) + aV;(t) + e(t) e(t) — y?a(t)* d(t) < 0 


2161 


(23) 


Equation (23) can be given as Eq.(24) by regarding Eq.(18) and Eq.(22), which may be represented in matrix 


inequality form as Eq.(25). 


e(t)! P; Baad (t) + e(t)! CT C;e(t)- 


e(t)" (A3 P; + Pj Aq; + aP,)e(t) +2 
) + d(t)* Dj, Daid(t) — y’d(t)"d(t) < 0 


2e(t)7 Ci Daid(t 


| nere ae S] ue CT Da 
0 


* * DI Da —?I | de 


One may obtain Eq.(26) by some simplifications in Eq.(25). 


I AN aP; Fi I 0 y 0 CF -7I 0 0 I e 
0 Bi P 0 Aci Belai I D} 0 I Ci Da 


By assuming Eq.(26) as NU ZNv, we have[32]: 


cli * 0 0 -| Aai Baai | <0 
0 BT j I cli cldi 
cldi * xo I+ DT Da; 0 I 


Thus, it can be concluded that: 


aP;+C?FC; P; C7 Dai 
* 0 0 <0 
* * —y! + D} Dai 


(24) 


(25) 


(26) 


(27) 


(28) 


Therefore, U may be achieved as [Acu — I Baai] taking into account that Ny columns are basis for the null space 


of U. Furthermore, according to [32], Ny and V are assumed as: 
I 0 
Ny-2|00!V«-|0I 0] 
0 I 


By using abovementioned results, Eq.(15) in Lemma 3 may be obtained as: 


Ai 0 aP,+CTC, P; CT Dai 
=I | X[OATO] + | AI | Xi[Aai — I Baai] + * 0 0 «0 (29) 
Bai 0 x * —y!I- DI Dai 


which can be cast into one LMI as Eq.(30). 


ISTE O12 O13 
LMIK(1) H * O22 O23 <0 
* x O33 


where 


O12 = Pi- XT + AAT Xi - ACT NF 
O15 = CT Dai + XP Bai — N; Dai 
O22 = —AX,; — AXP 
055 = AXP Bai — AN; Dai 
035 = DT Da; — ol 
o=% Ni = X] Li 


(30) 
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32. Fault sensitivity performance in the output estimation error of the switched observer 


By considering the Lyapunov function as Eq.(31) for faulty case, its derivative may be obtained as Eq.(32). 


V(t) = e(t)” Prie(t) 


V(t) = et)? Prie(t) + e(t) Prie(t) = elt)" (Ao; Pri + PriAai)e(t) + 2e(t) PriBasif (t) 


According to Remark 2, we need e, (1)? e, (t) for using lemma 1, which can be given as Eq.(33): 


e(t). e(t) = elt)" C; Cie(t) + 2e(0) C? Dg f(t) + IO DED pif) 
According to Remark 2 and Lemma2, one may obtain that: 
e(t)! (AT; Py; + Pri Acii)e(t) + 2e(t)" Pri Bayi f(t) + ae(t)? Prie(t)— 
e(t)" C? Cie(t) — 2e(t) C? Dyi f(t) — SOT DED if (t) + B f(t). f(t) «0 
which can be arranged into Eq.(35). 
e(t)? (A3; Pri + Pi Aai + oPyi)e(t) + 2e(t)” PriBagif(t) — e(t)! C7 Cie(t)— 
2e(0 C? Daf) + f(t)" (Df; Dri + PDS <0 


The matrix inequality representation of Eq.(35) is given as Eq.(36). 


AT Pri + PriAcii + aP ei PriBai 4 -CIC, -C7 Dji uj 
* 0 * BI -— D, Dji 
Equation Eq.(36) is given as Eq.(37) in form of NU ZNy. 
I AT. 0 a Pri - CICi Pri -C7 Dji I 0 
| 0 Bt —I | 0 0 0 Aci Bafi « 0 
el ft * 0 @I-— DF, Dji 0 —I 
Thus, Z matrix may be obtained as follows. 
a Pri = QUO Pri -CTIDg 
* * ggl- DD si 
U = [Acai —I Basil, V = [0 AI 0] 
Using Eq.(15) in Lemma 3, we have: 
AT; 0 &Pr TR GIG, Pri -CIDg 
—I Xj (0 AI 0] +| AL XiļAcu —I Bai] + * 0 0 
Bari 0 * * BI- D3} Dfi 


Altogether, a linear matrix inequality as Eq.(40) is achieved by assuming 7 = 37. 


V V1» V5 
LMI(2) i * V2» Wo3 «0 
* * YV353 


where 


Via = Pri — XT AAT X; — ACT NE 
Vi = -C7 Da +X] Bp -— N;Dji 
Woo = —AX; — AXP 
Vo; = AXT Br; — AN; Dfi 
V3 = TI — D4 Dji 


in which, 7, N;, Pr; and X; are unknown variables of the problem. 


«0 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 
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3.3. Simultaneous Optimal Fault sensitivity and Disturbance Attenuation Level 


According to Eq.(30) and Eq.(40), the disturbance attenuation level and fault sensitivity are defined in two 
separate LMI problems, which must be considered simultaneously in the observer design at a given disturbance 
attenuation level or fault sensitivity. In this paper, by defining a weighted LMI optimization problem, the optimal 
values of the disturbance attenuation and fault sensitivity are obtained. The theorem is given as follows. 

Theorem: Consider the continuous time switched system as Eq.(1)-Eq.(2) with ADT along with the 
switched observer as Eq.(5). The state estimation error dynamic Eq.(7) is asymptotically stable and satisfies the 
performance criteria as Eq.(9) and Eq.(10) for any nonzero d(t) € l2[0, oo), if there be the scalars values 7 > 0, 
u > 0 and the matrices P; > 0 Q; > 0 X; Nj such that the subsequent weighted LMI optimization problem has 
solution. 

m 


Us be TER LES (1— w)r 


s.L.DMI(1) «0, LMI(2) < 0 
P; < Pj, Pri < LP; 
where 0 € w < 1 is the weighting factor for simultaneous fault sensitivity and disturbance attenuation perfor- 
mances and its value is defined by minimization of the following cost function: 
B 


Lastly, the gains of the observer for different modes of the switched system, optimal disturbance attenuation level 
and optimal fault sensitivity of the output estimation error of the observer are achieved as: 


Li =X; 1 Ni, y= Va, B — VT 


3.4. Global Analytical Redundancy Relations in Error Form Using BG method 
The generic form of a GARR may be given in Eq.(41) [34]. 


GARR : v(y(9,...,,u9,...,u,0,a) (41) 


where y, u and 0 are denoted for outputs, inputs and parameters of the system, respectively. As in Eq.(3)-Eq.(4), 
a = [a4 a2... az] is a binary vector that determine the active mode of the switched system.. 

The GARRs are differential equations, which their orders are dependent on the sensor locations in the 
system. The lower order GARRs are more efficient for monitoring of a system. The generic form of the GARRs 
may be presented in matrix form by maximum second order assumption of derivatives. 


GARR : (9j, ý, y, à, à, u,0,a) = Mi(8, a)j(t) + M3(0, a)j(t) + Ms(0, a)y(t)+ 


Z,(6,a)ü(t) + Za(0,a)ü(t) + Z5(0,a)u(t) P 


The effects of disturbances, noises and parametric uncertainties on the GARRs may lead to false or missed 
alarms in the FD system. In this paper, a new method by combination of the BG method and observer method is 
proposed to overcome these problems in the FD system based on the BG method, which further leads to some 
benefits in the FD system. 

In normal case i.e. in the case that there are no disturbances, faults, noises and uncertainties in the system, 
the estimated states and outputs of the observer converge to the states and outputs of the system. Therefore, the 
GARRS by using estimated outputs may be given as: 


GARRDA : (ij, 9, i, ti, à,u,0,a) = Mi(0,a)j(t) + Ma(0, a)g(t) + Ma(0, a)g(t)+ 
Zi(0,a)ü(t) + Z2(0, a)u(t) + Z3(0, a)u(t) 


By subtracting Eq.(43) from Eq.(42), the GARRs in the error form, which is called EGARRs, are obtained. 


(43) 


EGARR : M,(0, a)éy(t) + Mo2(0, a)é,(t) + Ma(0, a)e, (t) (44) 


The output estimation errors in Eq.(44) come from the observer, which is designed in such a way that be 
more sensitive to faults and simultaneously the effect of disturbance is attenuated. Therefore, the effects of the 
disturbances are reduced in the GARRs, while the fault sensitivity is enhanced. 
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Flow 
source 


Tankl Tank2 


Figure 1. The two-tank system 


3.5. Robustness of the FD system against parametric uncertainties 


The fault detection and isolation are based on the EGARRs, which are dependent on the parameters of 
the BG model. In order to obtain a robust FD system against parametric uncertainties, the proposed method in 
[23] is considered in this paper. In robust fault detection based on the BG method, a GARR can be decoupled into 
nominal (GARRn) and uncertain part (UGARR) as follows. 


GARR = GARR,(a,0,,u,y) + UGARR(a, 69, u, y) (45) 


where 0, is the nominal part and dg is relative deviation compared to the nominal value of the parameter 0. 

The nominal part is the obtained GARRs, which contain the nominal values of the parameters and made 
the residuals. The uncertain part contains the uncertain parameters of the nominal GARRs, which further deter- 
mines the upper and lower bounds of the residuals. 

An EGARR may be given as Eq.(46) regarding this fact that the EGARRs are not directly dependent on 
the inputs and their distribution matrices. 


EGARR = EGARR,(a, 0n, ey) + UEGARR(a, 69, ey) (46) 


The adaptive thresholds on the residuals are defined as UEGARRs, which are dependent on the parametric 
uncertainty (ôg), mode of the system (a) and the output estimation error of the observer (ey). 

As mentioned above, the proposed method is not directly dependent on the inputs and distribution matrices 
of the inputs (Z4 to Z3). Thus, the robustness of the proposed method may be improved. 

Finally, the residual evaluation function and thresholds for fault detection are assumed as Eq.(47) and 
Eq.(48), respectively. 


lo+L 
In (t) = |Irell2 = (X rere)? (47) 
lo 
Jin = sup Ii (t) + |lexlla (48) 


d(t),u(t)€la,f (t)=0 


in which Jz (t) is the residual evaluation function, J;;, is the threshold on the residual evaluation function, lo is the 
initial time, L is the window length for residual evaluation function calculations and ex is the effect of parametric 
uncertainties on the residual, which is obtained from UEGARR. Finally, the faulty or fault free conditions are 
declared by comparison between the residual evaluation function and considered threshold. 


4. RESULTS AND DISCUSSION 

Consider the two-tank system in Fig.1, which involve a flow source, two output valves named as R; and 
FR» and one interconnection valve (R3). The valves could be in open or closed state, which determined the active 
dynamic of the switched system. Three configurations as Table I are assumed [35]. 

The BG models of the understudy two-tank system in integral and derivative causality in 20-sim software 
are shown as Fig.2 and Fig.3, respectively. 

As mentioned earlier, the valves could be in open or closed state. The coefficients a1, a» and ag are used 
to show the valve states of R1, Rg and R3, respectively, which are set to "0" in closed state and "1" in open state. 
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Table 1. Considered modes of the two-tank system 


2165 


Mode Ry Ro Rs A= [ay a2 às] 
Model off on off [0 1 0] 
Mode2 on on off [1 1 0] 
Mode3 off on on (0 1 1] 
c:cTl |B R:R2 C:CT2 B 
2 $ a2.f7|; 9 ^ 
Sent "Y "M. an tt 0 
3 alfa "as ft 
3t x 
tc 103 
*lal.f4 7la3.fD 
P t 
R:RI R: R3 


Figure 2. BG model of two-tank system in integral causality. 


In hydraulic domain, the effort and flow variables are regarded as pressure and flow, respectively. The subscript 
c is considered to determine the controlled junctions of the model. Three controlled junctions are assumed for 
three valves. In BG theory of hybrid systems, the flow variable of the connected bonds to the i^^ 1-type controlled 
junction is multiplied by a;. Two pressure sensors are used to measure pressures of the two tanks which is given 
as P, and P5. 

The state space representation of a hybrid system may be obtained in compact form using the BG model. 
The states of the system are considered as generalized momentum of inertia elements (pr) and generalized dis- 
placement of capacitive elements (qc). Hence, the state vector of the two-tank system is achieved as [qo qo] 


The flow through the valves is assumed laminar. This assumption leads to linear equations of the system 
[35]. Using the state space derivation method in [34], the governing equations in compact form are: 


ay a2 a2 
| + din (49) 
RıCrı Cn” nen 1 
+f ag 3 
J 


RəCr2 R3Cr3 
Yi = K 92; Y2 = 


- 99 
Cri Cra 


The parameters definition and their values are given in Table II. 

It is also noted that the tank capacity is defined as Cr = A/g wherein g = 9.81m/s? is the acceleration 
of gravity. 

Thus, the system matrices for three considered modes of the two-tank system are: 


go = ( 


F a2 ) 
q9 = q2 qo 
R2Cr1 


A= | —0.0007 0.0011 | Ge | 0.1962 0 | 


0.0007  —0.0011 0 0.3270 
à De:P2 
ccm DePl gg — ccm 7 
T T 
L7 22.07. s 7 
A rs 1 6 8 
St Qia n a2.f6 "Ut a2.f8 j 
3|al.f "| a3.f11 
£ £ 
1c 103 
al.f4 2) a3.f12 
il 
R:RI R:R3 


Figure 3. BG model of two-tank system in derivative causality. 
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Table 2. Parameters of two-tank system 


Variable Definition Value 
Aj Cross sectional area of tank 1 50m? 
Ag Cross sectional area of tank 2 30m? 
fü Output pipe of tank 1 resistance 300s/m? 
Ro Interconnection pipe resistance  3005s/m? 


Rs Output pipe of tank 1 resistance 100s/m? 


4 0.6 | 04 0. 
weight weight 


Figure 4. Disturbance attenuation level and fault sensitivity versus weighting factor 


Ay =| -0.0013 0.0011 | œ, _ [ 0.1962 0 
27 | «0007 -—0:00f1. | 94—- 0 0.3270 
A, =| -0-0007 0.0011 | o [0.1962 0 
3 | 0.0007 -—0.0044 |^ ^3 — 0 0.3270 


1 
By=By=Ba=| 4 | 


The other matrices of the system as in Eq.(1) and Eq.(2) are assumed as follows. 


1 0 
Ba = Ba = Bas = B1, Bp = Bf2 = Bf3 = | | 


0 1 
T 03 0 
Da — Da - Da - | 91 |. Dn 7 P - Dri - | 0 ea 


The GARRSs of the system may be obtained using the Causality Inversion Method (CIM) in [34]. 


dP, ag a1 
GARR, =C P, — Po) 4 P inc 50 
1 nU RC 2) [E q (50) 
dP» ag a3 
GARR = Č P — P. P,-—0 
2 T2 Re 1 5) + Rs 2 
The matrices M; and Z; can be calculated as: 
C 0 
M =0, Ma = | a es [t= ho 
T | m T 1 
My =| iene a 26 IE | 
m. Mats 0 


The upper bound of parameter uncertainty is assumed as 2% of their nominal values [29]. 
Solving the weighted LMI optimization problem for a = 0.08, u = 1.2 and \ = 0.9, the disturbance 


attenuation level and fault sensitivity versus weighting factor are given in Fig.4. 


By minimization the cost function of the Theorem, we have: 


w = 0.8 
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Figure 5. Switching signal 


The disturbance attenuation level, fault sensitivity and ADT are also given as: 
B = 0.4116, y = 0.2941, 7; = In(u)/a = 2.279 
One may obtain the following observer gains for three defined modes of the switched system. 


ie 2.0887 1.3471 | Eom | 2.0880 1.3512 | ds | 2.0789 1.3987 


—0.0522 0.3873 —0.521 0.3864 —0.0453 0.3514 
The disturbance signal in the switched system is considered as follows [32]. 
d(t) = (0.07 + 0.5e~°* cos 0.77t) u(t) (51) 


The switching signal is plotted in Fig.5. The EGARRs may be calculated as: 


de 1 2 ay 
EGARR; = Cr a tm (ey1 — €y2) 4 men 0 (52) 
déy2 az a 


3 
EGARR» = Cro di Rz (ey1 ey2) + R; €y2 = 0 


Thus, the uncertain parts of the EGARRs, which are considered in the thresholds on the residuals, are 
obtained as: 


de a a 
& = [don Cri qr |+ ba p (en 7 ew)l + Ba p enl (53) 
déy2 ag a3 
& = [fen Cra— l 4 i m, oul evi + lon. Rl 


The upper bound of thresholds due to the parametric uncertainties is defined by evaluating ||€;||2 for each 
residual. 

The simulation results are obtained for the case that the uncertainties in the parameters (maximum 2 96) 
are included in the system, but the observer and EGARRS are based on the nominal values of the parameters. At 
first, the first part of the thresholds according to disturbances is obtained in fault free case. The values for the 
thresholds are obtained as 0.0434 and 0.0063 for rı and ro respectively. Second, the thresholds on the residual 
evaluation function are achieved by calculating the parametric uncertainty parts i.e. ||e;.|]o . 

The obtained residuals and the residual evaluation function for fault free case in the presence of distur- 
bance and parametric uncertainties are shown in Fig.6 and Fig.7, respectively. 

As can be seen from the figures, the residual evaluation functions are not exceeded from their two-part 
thresholds (fixed and adaptive) and thus, no fault is declared. 

A fault signal (f1) is considered as a pulse signal with unit amplitude that occurred in 40 seconds and 
remains until 50 seconds for faulty situation. The residuals and residual evaluation functions are plotted in Fig.8 
and Fig.9, respectively. It can be observed from the figures that the fault sensitivity is enhanced in an efficient 
manner and can be efficiently discriminated from disturbance. The residual evaluation function of the first residual 
is exceeded from the defined threshold, which declares a faulty situation in the system. It is also clear that the first 
residual is sensitive to f; fault. 
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Figure 6. Residuals in fault free case 
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Figure 7. Residual evaluation function in fault free case 
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Figure 8. Residuals in faulty case 
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Figure 9. Residual evaluation function in faulty case 
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Figure 10. Residual evaluation function in conventional method for fault free case in the presence of parametric 
uncertainties 


In comparison to the case that the parametric uncertainties are not considered in the FD system, the 
proposed method leads to robustness of the FD system. The method without two-term threshold may lead to 
false alarm, which is shown in what follows. The obtained results for the non-faulty case in the presence of 
disturbance and parametric uncertainties are shown in Fig.10. The considered thresholds are based on the effects 
of disturbances in the system and the adaptive terms due to parametric uncertainties are not included. 

As can be seen from the figure, the residual are exceeded from the considered thresholds and occurrence 
of a fault is declared, while there is no fault in the system. The residuals for the proposed method may be given in 
Fig.9 in which there is no exceeding from the thresholds in this case. 

According to the simulation results, the proposed FD system is robust against parametric uncertainties 
and simultaneously attenuates disturbance level and enhances fault sensitivity. These criteria are achieved using a 
two-stage FD system by combination of a robust switching observer and BG method. 

The effects of the parametric uncertainties may lead to false alarm in the fault diagnosis system, as their 
effects are not included in the thresholds. The result in this case is shown in Fig.10, which declares a fault 
occurrence in comparison to the proposed method as in Fig.9. In fact, the new defined term as the effect of 
disturbance on the residuals will improve the robustness aspect of the FD system. The effects of disturbances are 
attenuated in the residuals by the switching observer in the first stage, which leads to smaller thresholds in the 
second stage. Smaller thresholds will also improve the efficiency of the FD system for fault detection in the case 
that the fault has small amplitude and small effect on the residuals consequently. 


5. CONCLUSION 

In this work, a new robust fault diagnosis system for continuous time switching system with average dwell 
time has been proposed based on combination of a switching observer and hybrid bond graph method, in which 
different criteria including fault sensitivity, disturbance attenuation and robustness against parametric uncertainties 
has been considered. Different criteria have been included in each stage of the two-stage proposed fault diagnosis 
system. In the first stage, a switching observer has been proposed that simultaneously enhances fault sensitivity 
and attenuates disturbance in the output estimation error of the observer in an optimal manner by defining and 
solving a weighted LMI optimization problem. Afterwards, the output estimation error of the observer has been 
used in a new form of the GARRs, which called EGARRs and are according to the output estimation error of the 
observer. The thresholds on the residuals contain two different parts including a fixed threshold due to the effect 
of disturbance and an adaptive threshold considering the parametric uncertainties of the systems in the EGARRs. 
A two-tank system has been considered and the simulation results have been shown the efficiency of the proposed 
method. 
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